Suppression of the critical current of a balanced SQUID 
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We present an experimental study of the magnetic flux dependence of the critical current of 
a balanced SQUID with three Josephson junctions in parallel. Unlike for ordinary dc SQUIDs, 
the suppression of the critical current does not depend on the exact parameters of the Josephson 
junctions. The suppression is essentially limited only by the inductances of the SQUID loops. We 
demonstrate a critical current suppression ratio of higher than 300 in a balanced SQUID with a 
maximum critical current 30 nA. 

PACS numbers: 85.25.Dq, 73.23.-b 
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Dc SQUIDs are routinely used to provide tunable crit- 
ical current Ic, e.g., in quantum computing applica- 
tions ^1, 2]. For example, some charge qubits and charge 
pumps would benefit if Ic could be tuned very close to 
zero 0, 0, Q. These devices require Coulomb blockade, 
and hence very small Josephson junctions (JJ) must be 
used. The range of the critical current of a dc SQUID is 
|/ci ± Ic2\ where Id are the critical currents of the indi- 
vidual JJs. It is not possible to fabricate two identical 
junctions, hence some residual critical current always ex- 
ists. As suggested in Refs. 0, 0, this problem can be 
completely eliminated by using a balanced SQUID, i.e. a 
structure with three parallel JJs in two superconducting 
loops with individual magnetic flux controls. Surpris- 
ingly only few experiments have utilized individual on- 
chip flux controls instead of a homogeneous external flux 
(see e.g. Refs. d, 0|)- In this paper, we present to our 
knowledge the first experimental study of the balanced 
SQUID. Our prime motivation to search for high critical 
current suppression ratio is to improve the accuracy of 
the Cooper pair sluice [3| . It is a current pump that can 
produce a current as high as 1 nA 0, Q , but whose accu- 
racy is not yet sufficient for a quantum current standard. 

A schematic picture of the balanced SQUID is pre- 
sented in Fig. [T^. It or any other system of n JJs in 
parallel forming n — 1 loops can be modeled as follows. 
The total current is the sum of the currents through 
the individual JJs, and the vector of the JJ currents is 
I = Ic sin 4>, where Ic is a diagonal matrix of the critical 
currents Id of the JJs. A vector containing the sines of 
the phases (pi over the JJs is denoted by sin 0. The effect 
of the magnetic flux $i on loop 1 is 
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Lii/ci sin^i + Li2/c2 sin(?!)2), (1) 



where $o is the flux quantum. The magnetic flux induced 
to the loop by the currents flowing in the branches j of 
the structure are proportional to the inductances Lij. 
For several loops, Eq. ([1]) can be written in a matrix 
form: M</> = f - LI. Here, f = (27r/<i>o)*, where * 
is a vector of the fluxes of the loops. The {n — 1) x n 
inductance matrix L contains the factor 27r/$o- 




FIG. 1: (a) Balanced SQUID. The middle junction is larger 
than the others, (b) Large detector junction in parallel with 
the balanced SQUID, (c) Scanning electron micrograph of 
the sample showing the on-chip coils and the narrow SQUID 
loops, and a simplified sketch of the measurement setup. Here 
the resistance in series with the pump is -Rbias = 100 kf2. (d) 
Magnified view of the junctions and a sketch of the SQUID 
loops. 



In the limit of zero inductances, the total current is ob- 
tained in an analytic form / = I^i sin((/)i) + /c2 sin((/)i — 
fi) + ... + Icn sin(0i - /i - ... - fn-i), where f., = 
27r(f>i/$o. The critical current of the whole structure 
/c(/i, . . . , fn-i) as a function of the magnetic fluxes can 
be solved with the help of trigonometric identities. The 
full model including the inductances requires extensive 
computation over the phases of the JJs, whereas the an- 
alytic model allows direct calculation of the critical cur- 
rent. Therefore, the analytic zero-inductance model is 
helpful for the analysis and parameter fitting of the mea- 
surement results when the inductances are small. 

Since we are interested in the suppression of a critical 
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current that is originally of the order of tens of nanoam- 
peres, we must be able to measure the critical current at 
least at 0.1 nA level. So low currents belong to the regime 
of phase diffusion, where a current-voltage {IV) mea- 
surement gives a maximum supercurrent roughly propor- 
tional to I^. Moreover, the result is very sensitive to 
the electromagnetic environment Q. Even an IV mea- 
surement of a critical current of about 50 nA requires a 
special environment Therefore we add a detector 

junction with 1^4 > 100 nA in parallel with our balanced 
SQUID, see Fig. [TJd. This superconducting shunt pro- 
tects the balanced SQUID from the environment. The 
critical current of the four-junction system exhibits sinu- 
soidal modulation around /c4 as a function of the flux 
$3. The amplitude of the modulation equals the critical 
current of the balanced SQUID. 

We measure switching into the normal state to de- 
termine the critical current at a certain magnetic flux 
set $ by feeding current pulses of constant length and 
variable height through the device. Switching produces 
a voltage pulse. We measure the switching probabil- 
ity at typically 5 points and determine the current, /50, 
at which the system has a 50 % probability to switch 
to the normal state. At the measurement temperature 
125 mK, with /c4 « 200 nA and with about 1 ms pulse 
lengths, /so = ale, where the factor a is between 0.6 
and 0.7 [ll|. To distinguish between the real critical cur- 
rents and the measured values obtained from the switch- 
ing experiments, we use the notation id = aid for the 
measured critical currents. Here the index i refers to an 
individual JJ or to the balanced SQUID. 

The measured flux modulation of the critical current is 
presented in Fig. [2^. For each data point of the 2D chart, 
we measured /50 for 11 different currents in coil 3 map- 
ping about one flux quantum in $3. The critical current 
ich of the balanced SQUID was extracted from the dif- 
ference of the maximum and minimum critical current of 
the whole structure as a function of the current in coil 3. 
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FIG. 2; (a) Measured critical current ich of the balanced 
SQUID as a function of the coil currents. The maximum is 
shifted from zero current due to an offset flux, (b) Respective 
theoretical flux modulation of the critical current calculated 
with the parameters fitted from the measurement. 

The critical currents of junctions 1 and 3 are almost 
equal in our balanced SQUID design, see Fig.[TJ Junction 



2 is larger. The critical current of the SQUID is tunable 
to zero if Ic2 < ^ci + Ic3- The flux modulation of the 
SQUID is periodic in the square with $i,<I>2 G [0,$o]j 
see Fig. [51 The critical current has the maximum value 
at the corners of the square. In the middle region, there 
are two minima. If Ic2 > Ici + /c3, there would be only 
one minimum with value Ic2 — [Id + Ics)- The two min- 
ima are close to the line (0, 0) ($0i ^o)- In the Cooper 
pair sluice, one could use a coil coupled symmetrically 
to the loops instead of the individual couplings exploited 
here. It would then be possible to move along the di- 
rection (0, 0) — > ($0j ^o) from the minimum close to the 
maximum with a single rf control. An asymmetrically 
coupled dc controlled coil would also be required, but dc 
signals are easier to implement. 

We fitted the zero-inductance model to the entire 3D 
flux modulation data. As fitting parameters, we used 
the critical currents of the junctions, ic, the offset flux 
vector ^offset; and the matrix K depicting the couplings 
and cross couplings between the coils and the SQUID 
loops: $ = ^offset + KIcoii. The fitted parameters were 
used to calculate the theoretical flux modulation chart 
presented in Fig. ^p. The resulting critical currents are 
id = 7.4 nA, ic2 = 13.5 nA, ic3 = 9.0 nA, and Zc4 = 
142.5 nA. The maximum critical currents of the balanced 
SQUID and the whole structure are about 30 nA and 
172 nA, respectively. 

Next, we performed a measurement at 17 x 17 points 
near the minimum of the critical current of the balanced 
SQUID with about $o/400 step in $i and $2- Since the 
coupling of the detector coil 3 to the balanced SQUID 
loops is expected to have the most pronounced effect near 
the minimum, we fitted the terms Ki^, and K23 of the ffux 
coupling matrix again to these data. The ffux coupling 
matrix combined from the two fits is 



K 



2.52 0.09 0.10 
0.04 2.39 0.11 
0.05 0.06 2.63 



<I>o/mA. 



(2) 



The cross coupling terms K13 and K23 are somewhat 
larger than expected, partly because the bonding wires 
of coil 3 pass near the sample. This enhances the effect 
that the fluxes $1 and ^'2 do not remain constant when 
the flux sweep over $3 is being performed. 

One of the measured flux modulation curves with re- 
spect to the detector coil 3 is shown in Fig. [3] by open 
circles. The chosen curve corresponds to the measure- 
ment point ($1, $2) that is closest to the fitted minimum 
of the critical current of the balanced SQUID. Due to 
the cross coupling, the lowest critical currents cannot be 
directly extracted from the modulation amplitude. The 
zero-inductance model gives the modulation curve shown 
by the black line and a critical current 0.035 nA at this 
measurement point. 

The parameter fits were performed to data extending 
over one ffux quantum in the detector loop. To test 
the validity of our model, we extrapolated the ffts over 
four periods in $3 and compared to the measurements. 
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assuming the resistivity of the aluminum loops to be 
p K, 2.1 /iilcm. The total calculated inductance matrix 
is 
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FIG. 3: Detector flux modulations (<l>i^PP'''^'^ — Kzil^oiVi) near 
the minimum of icb • The horizontal lines present change from 
ic4 as a percentage of the maximum of icb- The open circles 
are an example of the data used for the parameter fits of 
the cross coupling terms and -K'23. The black line is the 
corresponding curve of the zero-inductance model yielding a 
critical current i^b = 0.035 nA. The data spanning four flux 
quanta presented by the red squares and the blue circles were 
not used for fitting. The corresponding red and blue lines are 
the predicted flux modulations yielding icb = 0.17 nA and 
icb = 1.8 nA, respectively. The inset shows typical results 
of similar critical current measurements of two ordinary dc 
SQUIDs, expressed as a function of the applied magnetic flux 
of the SQUID. The minima of the critical currents are 2.5 % 
and 4.6 % of the maxima, which are about 30 nA for both 
SQUIDs. 



The theoretical predictions and the measured data at two 
((f>i,$2) points near the minimum are shown in Fig. [31 
The only fitting parameter was the overall magnitude of 
each data set. The deviations of the measurement points 
from the theoretical curves are less than about 0.2 nA, 
and no systematics of the deviations is observed. 

The remaining sources of residual critical current in 
the balanced SQUID are noise and the inductances of 
the loops. The effect of the inductance L is proportional 
to the magnetic flux that a current circulating in the loop 
can produce: 27rL/ci/$o, see Eq. ([1]). This effect can be 
dominating with large JJs, which can be fabricated with 
smaller relative parameter scatter. Hence the balanced 
SQUID does not outperform the dc SQUID in that case. 

The inductances of the SQUID loops consist of geo- 
metric and kinetic inductance. The geometric induc- 
tances were calculated with the superconducting ver- 
sion of FastHenry fl^, and the kinetic inductance by 



L/(27r/$o) = 









pH. (3) 



Here we have neglected the cross coupling terms except 
^21, which is large due to the physical position of current 
injection to the structure, as shown in Fig. [TJ 

We simulated the flux modulations of the critical cur- 
rent with inductances from L of Eq. (3) up to 20L, and 
with or without the detector junction. The critical cur- 
rent suppression ratio of the balanced SQUID is about 
500 with L and about 350 with 2L. In the four-junction 
model, increasing the inductance to 5L or higher breaks 
the simple sinusoidal modulation, and the modulation 
amplitude is increased. Such effects were not observed in 
our measurements, which supports the estimated values 
of the inductances. 

Noise is expected to round off the apparent critical cur- 
rent suppression. An obvious noise source is the global 
magnetic field, which has been diminished by a high-/i 
and a superconducting shield. In low-Tc SQUIDs, the 
main 1// noise sources are the critical-current fluctua- 
tions and magnetic flux noise related to defects close to 
the superconducting loop lines [31 . Our measurement 
is, however, relatively insensitive to magnetic flux noise, 
since the phase of the detector flux modulation turns 
by 27r around the minimum of the critical current of the 
balanced SQUID. Thus, magnetic flux noise creates mod- 
ulation signals with all possible phases, and in the first 
order, this modulation caused by the noise averages to 
zero. The effect of flux or critical-current noise was not 
observed in the measurements. 

In conclusion, we have demonstrated experimentally, 
that the critical current of the balanced SQUID agrees 
well with the zero-inductance model. This is supported 
by theoretical estimation of the effect of the inductances. 
A conservative estimate for the critical current suppres- 
sion factor is 300 or more. This number is an order of 
magnitude higher than the typical values obtained for 
two-junction dc SQUIDs with similar maximum critical 
currents. Furthermore, the suppression ratio does not 
essentially depend on the exact critical currents of the 
individual junctions, and hence the fabrication yield is 
high. 
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